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We consider the problem ol tunneling between spinless ID Luttinger liquid chains. We show how to 
map the problem onto the 4-state chiral clock model together with a free boson for the total charge 
mode. We use a variety of results, some of them exact, from the integrable chiral Potts model and 
the study of commensurate/incommensurate transitions to deduce the physics of coupled Luttinger 
liquids. For those intrachain interaction strengths for which interchain tunneling is relevant, we 
find that it can lead to the formation of symmetric and antisymmetric bands with split Fermi 
surfaces, depending on the relative strengths of interchain tunneling and interchain interactions. 
With split Fermi surfaces, interchain transport is coherent. It is not possible to have two gapless 
Fermi surfaces with the same Fermi momentum when interchain tunneling is relevant. However, 
interchain interactions can drive the formation of a gap, in which case interchain transport will 
be incoherent. We comment on the possible relevance of our results to c-axis transport in high-T c 
superconductors. 

PACS numbers: 72.27.+a, 72.10Bg,74.25.-q 
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I. INTRODUCTION 

The problem of tunneling between II? Luttinger liq- 
uids has been a focus of interest and controversy. Much 
of the interest stems from the unusual phenomenology of 
the cuprate superconductors. The normal state of these 
materials exhibits insulating (dp c /dT < 0) or, at best, 
very poor c-axis transport coexisting with metallic ab- 
plane transport IB. Moreover, angle-resolved photoemis- 
sion [^1^] reveals a Fermi surface with no dispersion along 
the c-axis, despite the fact that the hopping matrix ele- 
ments are quite substantial between the two planes of a 
Bi2212 bi-layer. Apparently, single-electron motion along 
the c-axis is very strongly suppressed while motion in the 
a6-plane is not. 

It has been suggested that this peculiar situation 
is characteristic of a putative non-Fermi liquid normal 
state. ID Luttinger liquids are the best-understood non- 
Fermi liquid metals; as a result, they are a natural test- 
ing ground for such ideas. Indeed, it has been suggested 
H that the type of behavior seen in the cuprates oc- 
curs when two ID Luttinger liquids are coupled by weak 
single-electron hopping between the ID chains. 

When intrachain interactions are extremely strong, in- 
terchain hopping between ID Luttinger liquids is irrel- 
evant in the renormalization group sense. Extremely 
strong means much stronger than can be reached in the 
infinite-?/ Hubbard model, and so strong that the fermion 
spectral function has no singularity whatsoever at the 
Fermi surface. In such a case, it is possible to have 
metallic behavior within a chain but insulating transport 
between chains. For weak and even moderately strong in- 
trachain interaction strengths, however, interchain hop- 



ping is a relevant perturbation. The natural assumption 
is that this leads to metallic interchain conduction. 

Many analyses P-[7|j9|,[To|,pt diagonalize the interchain 
hopping term first, forming symmetric and antisymmet- 
ric bands with displaced Fermi surfaces in the 2-chain 
case. Interactions are then treated as perturbations. Of 
course, if these perturbations could be treated exactly, 
the resulting solution would be the same as that which 
would be obtained by diagonalizing intrachain interac- 
tions first and then solving interchain hopping (which we 
will do non-perturbatively and, in one interesting limit, 
exactly, as we describe below). However, this need not 
be true for approximate solutions; since the approximate 
solutions of [p|-[7[P,^0|,p[| assume c-axis dispersion (i.e. 
dispersion in the direction perpendicular to the chains, 
which we'll call the c-axis) by taking the existence of sym- 
metric and antisymmetric bands as a starting point, it is 
not surprising that they do not find the scenario of Q . 

Consequently, it has not been possible to rule out 
claims Q that interchain hopping, though relevant, does 
not lead to coherent transport (we discuss the possible 
meaning of this below) and, thus, does not lead to c- 
axis dispersion or metallic interchain transport. These 
claims have been supported by appealing to the spin- 
boson problem - in which a two-level system is coupled 
to a bath of bosons - which has a parameter regime in 
which tunneling between the two levels is relevant, but 
there are no coherent oscillations ]l2| , p^| . The authors 
of Q have claimed that there is a similar range of intra- 
chain interaction strengths for which interchain hopping 
between ID Luttinger liquids is relevant but does not 
lead to metallic interchain transport. 

To resolve this issue, we must address whether there 
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is c-axis dispersion and metallic interchain transport. If 
the minimum energy for adding an electron occurs at two 
different Fermi surfaces, ±k F and ±k F , corresponding to 
the symmetric and antisymmetric bands (i.e. if there are 
two different wavcvcctors for the minimum energy exci- 
tations which are even and odd under the Z 2 symmetry 
which exchanges the two Luttinger liquids), then there 
is c-axis dispersion. If, in addition, this minimum en- 
ergy is zero, then interchain transport will be metallic. 
In this case, various long-wavelength power-law correla- 
tion functions will oscillate at wavevector ± (k% — k F ). 
According to there is no gap, but this splitting does 
not occur: k s F = k F . 

However, the problem is complicated somewhat by the 
fact that interchain interactions can drive the formation 
of a gap. When this occurs, it is possible to have c-axis 
dispersion without metallic interchain transport. In this 
case, the minimum energy for adding an electron occurs 
at two different 'Fermi surfaces', ±k F and ±k F , but this 
energy is non-zero. The opening of such a gap would 
not explain the absence of a bi-layer splitting in photoc- 
mission studies of the cuprates (|J^] since it occurs 'on 
top' of the bi-layer splitting. On the other hand, it is 
also possible for the gapped phase to show no c-axis dis- 
persion, so that the minimum energy (which is non-zero) 
for adding an electron occurs at a single 'Fermi surface'. 
Clearly, there is no phase transition between these two 
limits of the gapped phase, but they embody a qualita- 
tive distinction which is important for photoemission in 
the cuprates. The purpose of this paper is to determine 
which of these possibilities is realized. We do this by di- 
agonalizing intrachain interactions first and then treating 
interchain tunneling and interchain interactions by non- 
perturbative methods, some of which are exact. 

The main technical difficulty which has prevented the 
resolution of the problem of tunneling between Luttinger 
liquids is that tunneling is a chiral perturbation - i.e. the 
interchain tunneling term transforms non-trivially under 
the 0(2) rotation group of 1 + 1-dimensional Euclidean 
(in imaginary-time formalism) spacetime. Physically, the 
tunneling operator is chiral because it mixes right-moving 
fcrmions in one chain with right-moving fermions in the 
other, and similarly for left-moving fermions; it does not 
mix right- movers with left-movers. This stands in con- 
trast to the case of a non-chiral perturbation - which 
mixes right-movers with left-movers - for which the 1 + 1- 
dimensional sine-Gordon model is the archetype. As we 
know from this exactly soluble models, non-chiral pertur- 
bations, if relevant, usually lead to the formation of a gap. 
Chiral perturbations are poorly understood, however. At 
the free fermion point, interchain tunneling is a chiral 
perturbation which is soluble because it is quadratic in 
fermion operators; it leads to the formation of symmetric 
and antisymmetric bands with Fermi surfaces which are 
shifted from the unperturbed Fermi surface. However, 
the behavior at this trivial point may not be generic. Un- 



til now, it has not been possible to appeal to a non-trivial 
exactly soluble model in which the fate of a relevant chi- 
ral perturbation is understood. 

In this paper, we consider spinless ID Luttinger liq- 
uids coupled by tunneling. We believe that spin is an 
inessential complication for the issue at hand; the claims 
of [Q apply to spinful and spinless Luttinger liquids, as 
do possible analogies to the spin-boson problem. We map 
this model to a statistical mechanical model known as the 
generalized 4-state chiral clock model. The phase struc- 
ture of this model is fairly well understood |l4],[l5| • This 
chiral clock model has a disordered phase and both com- 
mensurate and incommensurate ordered phases. As we 
discuss below, the incommensurate phase corresponds to 
an electronic state with c-axis dispersion and metallic in- 
terchain transport. This picture receives further support 
from the consideration of a particular limit in which the 
model becomes the integrable chiral Potts model which, 
according to exact results fT^JT^| , is in an incommensu- 
rate phase. 

In the interesting regime in which interchain tunnel- 
ing is relevant, there are two possible phases. Depending 
on the relative strength of interchain interactions and in- 
terchain hopping, the system can either be in a gapless 
incommensurate phase or a gapped phase. In the incom- 
mensurate gapless phase, the coupled Luttinger chains 
exhibit c-axis dispersion and metallic interchain trans- 
port despite strong intrachain non-Fermi liquid proper- 
ties. It is not possible for the system to be gapless and 
commensurate; in other words, there is no gapless phase 
with two identical Fermi surfaces if the tunneling matrix 
element is non-vanishing. Hence, some of the claims of 
H are incorrect. The gapped phase, on the other hand, 
may or may not have displaced Fermi surfaces. Thus, it is 
possible that the mysterious c-axis physics of the high-T c 
cuprates can be explained by a simple analogy with the 
ID Luttinger liquid, but only in the somewhat extreme 
parameter regimes in which c-axis tunneling is irrelevant 
or a gap forms which is strong enough to preclude c-axis 
dispersion. 

In section 2, we introduce the model of coupled Lut- 
tinger liquids and give its bosonized formulation. In sec- 
tion 3, we review the results for tunneling between Fermi 
liquids, and see how incommensurability arises there. In 
section 4, we map the Luttinger problem to a generalized 
4-state chiral clock model. In section 5, we use results 
from the chiral clock model to understand the physics of 
coupled Luttinger liquids. 

II. THE MODEL 

The Luttinger liquid is a model of interacting gapless 
fermions in one spatial dimension. At the microscopic 
level, it might be described by the lattice Hamiltonian 
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where i, j label sites along the ID chains, Cj annihilates 
a fermion at site i, and rii — CiC^ In the long- wavelength 
limit, the corresponding Luttingcr liquid action is: 



where the fermion creation operator is 

c(x) = e~ ikpx c R (x) + e lkFX c L (x). 



(2) 



(3) 



As is well known (see |18| and references therein), the 
Luttinger model may be rewritten in terms of a single 
free boson with action 



Sri = dx dt 



i [{d t pf-v\d x pf) (4) 



The field, <p, is taken to be an angular variable satisfy- 
ing the periodicity condition ip = ip + 2ixr. Here r = 
^/g 7 where g is determined by the interaction strength, 
g = \J{v — u)/(v + u). The coupling g < 1 for repul- 
sive interactions, g > 1 for attractive interactions. Our 
results apply to both cases because the Hamiltonian is 
invariant under the duality symmetry g — > 1/g. In the 
long-wavelength limit, the theory is 'Lorentz invariant', 
with v playing the role of the speed of light. 

The operators of the theory are defined in terms of 
chiral fields obeying (dt ± vd x ) pr.l = 0. Then ip = pl + 
ipu, while the dual field (p = (pi, — ipR. It is convenient 
to label the operators by their "electric" charge m and 
"magnetic" charge n, defined as 



r i ( ■ m ■ ~\ 

{m, n\ = exp \ i—ip + inr<p j 



(5) 



The left and right scaling dimensions of this operator are 
then 
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This means that when r = ^/g, the overall scaling dimen- 
sion x = fiL + hfi and Lorentz spin s = — of the 
operator {m, n} are 



x = — h ?7 q . 

Ag y 



The electron operator in particular is the {1, 1/2} oper- 
ator, namely 



c L = e 



(7) 



Then c[ = {-1,-1/2}, while c R = {1,-1/2} and 
c R = { — 1,1/2}. The electron operators have Lorentz 



spin ±1/2, meaning that they have non-trivial proper- 
ties under Lorentz transformations (i.e. rotations in the 
in 1+1 dimensional spacetime). Operators with non-zero 
Lorentz spin are called chiral. Adding a chiral operator 
to the action breaks Lorentz invariance. 

In this paper, we consider the effects of coupling two 
identical spinless Luttinger liquids. There are two types 
of couplings which will concern us here, interchain hop- 
ping and interchain interactions. Interchain hopping 
takes the form: 



-t± 



h.c. 



while interchain interactions take the form: 
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i - J\)ni,iTij,2 



(8) 



(9) 



Here, c i 1 annihilates an electron at site i of chain 1 while 
c,j 2 annihilates an electron at site i of chain 2; 77^1 and 
77^2 have corresponding definitions. 

In the long- wavelength limit, we can rewrite the model 
in terms of two free bosons ip± and if2- However, the 
interactions depend only on the difference of the two, so 
we define the linear combinations 
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V2 



(ip! ± tp 2 ) 



(10) 



Then interchain hopping takes the form: 



£ c = -tf e 75[(v^+ 1 /\/ff)y--R-(\/9- 1 /v / 9)¥'-- L ] + c.c. 

-t\ e 7f[(v / 9+ 1 /v / 9)v-L-(v^-l/v / 9)¥>-n] + c( , /-q) 

In the decoupled theory, the dimensions just add, so this 
operator has dimension and spin twice as large as the 
fermion, namely x = (g + l/g)/2 and Lorentz spin ±1. 
Thus this theory has chiral operators in the action and 
breaks Lorentz invariance. In a model of coupled chains 
which is parity- invariant, t\ = t±. However, it will be 
convenient to allow these two hopping parameters to be 
independent so as to separate various effects. In partic- 
ular, the case t\ = 0, ^ will be of special interest. 
Interchain interactions take the form: 



£v ± = V 2kF cos(V2^-) + JcoB(y/2/g<fu). (12) 

The first is a pure interaction term of dimension 2g: no 
charge is transferred between the two liquids. The second 
is a Josephson tunneling or pair hopping term; it trans- 
fers a pair of electrons from one chain to the other. Even 
if we set these terms to zero by hand (and there is no 
reason to do so) , second-order perturbation theory in the 
interchain hopping induces them. This model is invari- 
ant under the electric-magnetic duality g — > 1/g: under 
this symmetry the pair hopping and interchain interac- 
tions are swapped. We could also add zero-momentum 
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density-density and current-current interactions of the 
form dx^ifii d x .t(p2, but they will not change the results 
described here: they only renormalize the coupling g and 
the Fermi velocity v, so we absorb them in redefinitions 
of the couplings. 

Thus the interactions of two coupled Luttingcr liquids 
can be described by a single boson </?_; the boson ip + is 
free and describes the total charge, which is unaffected by 
the interchain coupling. It is useful to rewrite the interac- 
tion terms in the action in terms of electric and magnetic 
charges with a radius r — y/g/2. The interchain hopping 
terms are of the form (Q), namely {±1, ±1}_, where the 
subscript indicates the extra \[2 in the relation of the ra- 
dius to g. The V^fej? term is {0,±2}_, while the J term 
is{±2,0}_. 

The lowest-order RG equations for these perturbations 
are: 



R,L 



O 



(13) 



^ = (2-2/g)J+0(tfti) 
^ = {2 -2g)V 2kp +0(t«tl) 

Several important features of the model are apparent in 
these equations. 

(1) For g ^ 1, either the pair- hopping or the 2fci? inter- 
action is relevant, and the other one is irrelevant. Thus 
when the electrons interact, there is always a relevant 
non-chiral term in the action. 

(2) For 1 > g > 2 - V3 = .268 . . ., tf, t L ± and V 2 k F are 
relevant. For g < 2 — Vo tunneling between Luttingcr 
liquids is irrelevant. For g < 2 — -\/3, the fermion op- 
erator has dimension 1, so its spectral function has no 
singularity at the Fermi surface (for 2 — \/3 < g < 1, 
there is at least a power- law divergence; for g = 1, i.e. 
free fermions, a S- function singularity); hence, this is the 
limit of extremely strong interactions. 

(3) If t\ = and V 2kF = J = but tf ^ 0, then all 
the higher-order terms in the RG equation for t± vanish: 
higher powers of t± are of higher Lorentz spin and thus 
cannot renormalize rj. This is in stark contrast to the 
pair-hopping and 2&f interaction terms which do renor- 
malize themselves. 

Our aim in this paper is to discuss the fate of the system 
under the flow of these relevant perturbations. 



III. WEAK-COUPLING PHYSICS 



We first consider the free fermion limit, g = 1. If 



we set V 2k = J = 0, then we simply have two Fermi 
liquid chains coupled by single- fermion hopping. This 
special case was discussed in detail in p9[. The hopping 



operators are of dimension (1,0) and (0,1). As a result 
of the interchain hopping, symmetric and antisymmetric 
bands are formed with Fermi surfaces displaced from fcj?. 
For the purpose of comparison with g ^ 1, we describe 
this in a little more detail. 
The action is 

S = dx dt( c\ RL (d t ± vd x ) c IRL + 

i=iX J ^ 

- if c\r c 2R + cc - t\ c\ L c 2L + c.c ) (14) 



where 7=1,2 specifies the chain. We diagonalize this by 
forming symmetric and antisymmetric bands in terms of 

c s , aR ,L( x ) = ( c i«,l(*) ± c 2R Jx)) e^" ^ (15) 
The action is simply: 

S = J dx J dt (c\ R L (d t ± v s d x ) c sRL + 

c Ir,l ( d t ± v a d x ) c aRJ )j . (16) 

If the fermion spectra are truly strictly linear for t^' L = 
0, then v s = v a = v. otherwise, these velocities will be 
shifted. 

It is instructive to re-cast these steps in the bosonic 
representation, which is the most useful representation 
for g ^ 1. Then, 

S = Jdxjdt-^- (W±) 2 -v 2 {d xV ±f) - 

dx J dt (tle 1 ^^ + t\e tv -^ + c.c.) (17) 

The SU (2) symmetry which interchanges the two un- 
coupled chains can be used to transform this into a more 
tractable form. Performing a rotation generated by the 
unitary operator 



U = exp (| (j*-*^ - , ' : " v - + L - /?)) (IS) 



?V2 _ g— 

we transform the action into: 



S= Jdxjdt^- ((9^'±) 2 - f 2 (cW±) 

dx / dt (t±d-(p'_ R + t^d+(p'_ L ) (19) 



where <p' is the rotated field, d± = dt ± vd x and, for sim- 
plicity, we have assumed that t^ L are real. The second 
line of ([n]) can be eliminated by making the change of 
variables: <p'_ R L — ► <p'-R l + ^J*x' ( v tT%)- As a result, 
various correlations functions of the unshifted variables 
will be spatially modulated: 



,V_R(aO/V5p-V-H(0)/V2\ 



(20) 
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The same is true of the unprimed fields. This spatial 
modulation of a correlation function is characteristic of 
an incommensurate phase, and can be taken as the def- 
inition of such a phase This behavior still occurs 
even if tf ^ t\. In particular, if tf ^ 0, but t\ = 0, 
the modulation occurs in the right-moving sector, but 
not the left-moving one. The relation to the fermionic 
representation is: 



(21) 



so comparing with (|) we see how in an incommensurate 
phase the fermi surfaces are shifted. 

To summarize, the hallmark of c-axis dispersion and 
metallic c-axis transport is the oscillation of the power- 
law correlation functions of (/?_, as exemplified by eq. 
(pp|). This is equivalent to the statement that gapless 
fermionic excitations occur at two different Fermi sur- 
faces. 

Let us now turn on weak interactions and move away 
from the free fermion point, but remaining in the weak- 
coupling regime. In this regime, it makes sense to con- 
sider both intra- and interchain interactions as a per- 
turbation of the system in which interchain hopping has 
already been solved (|l6|). These interactions take the 
form: 



dx / dt 



u s c Ir c sR c Il c sL + 
U a C aR c aR C aL C aL 



C Ir C sR. c \l C sL + c \r c ciR c Il C sL 



+ U X ( 

+ U t (cIr C Il C aR C a L + 4«4 



(22) 



The renormalization group equations for these interac- 
tions read: 

du KM , V., + V„, 9 



dt 
du x 

~dJ 

du t 

~dl 



2v x 



2v s 



2v a 



2u, 



u, 



(23) 



For v s = v a , these are the Kosterlitz-Thouless equa- 
tions for Ut and u s + u a — 2u x . In the limit of small 
Ut, the system is described by two gapless decoupled 
Luttinger liquids (corresponding to the symmetric and 
antisymmetric bands) if u s + u a — 2u x > 0. If, in the 
same limit, we instead have u s + u a — 2u x < 0, then a 
gap is opened. The action ( p2[ ) is Lorentz invariant in 
this limit, so this gap cannot affect the splitting between 
the Fermi momenta of the symmetric and antisymmetric 
bands. The lowest energy single-fermion excitations will 
be at « kp ± ^jtj- When we move away from v s = v a 
(which will certainly happen if the interaction strengths 
are different in the two bands), we can no longer appeal to 



Lorentz invariance, and it is possible that the symmetric 
and antisymmetric bands are not split in momentum. In 
the limit of very large u s + u a + 2u x , the entire approach 
breaks down because interchain tunneling between the 
two chains is irrelevant. In this limit, symmetric and an- 
tisymmetric bands do not form, and c-axis transport is 
insulating. In subsequent sections, we address the ques- 
tion of whether this can happen in the strong-coupling 
limit even when interchain tunneling is still relevant. 



IV. THE CHIRAL CLOCK MODEL 

In this section we introduce a lattice model which, in 
the continuum limit, is described by the above coupled 
Luttinger liquids, including the chiral operators. In the 
next section we will utilize known results on the incom- 
mensurate phases of this lattice model to extract physics 
for our problem. 

The lattice model is a generalized chiral clock model, 
with nearest-neighbor interactions on a square lattice. It 
is written in terms of discrete "clock" variables 2air/N, 
where a is an integer ranging from to N — 1. The 
interaction energy between neighbors of values a and b is 
required to have a Z„ symmetry and, thus, is of the form 
£ (a, 6) = £(a — b) — £(a — b + N). The general action is 



N-l 



cos 



2tt. 

— {jD x a-Aj) 



-Kj cos 



2tt. 
~N 



(jD y a - Aj 



(24) 



where D is the lattice derivative, namely D x a = a(x + 
1, y) ~ a(x, y) and D y a = a{x, y + 1) — a{x, y). We can 
obtain the physics we want by taking 



'v. 'v. 



K N -j = Kj 



Ai 



-A N _! = A A x = 
A, = A,- = 



-A N ^ = A 
otherwise 



We will also take Kj to be real, but A and A need not 
be real. Usually, clock models have Kj — for j > 1, 
but it will prove very useful to relax this condition. This 



model was originally introduced in |14|, which discussed 



the special case Ki = 1/T and Kj = for j > 1, and 
A = 0. When A =/= 0, the model is chiral. For < 
A < 1/2, it is energetically more favorable for D x a to 
be 1 than for it to be —1: a tends to rise (cyclically) as 
x moves to the right. The fact that at low temperature 
there are two possibilities (no rise or a rise of +1) leads 
to the existence of an incommensurate phase, to which 
we will return in the next section. 

For couplings close to a critical point, one expects the 
continuum physics of the chiral clock model to be de- 
scribed by a field theory. For A = A = 0, a number of 
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critical points are well understood. For N=3, there is 
a single critical point at K\ = [ln(\/3 + l)]/3, the S3- 
symmetric 3-state Potts model critical point. The cor- 
responding field theory is a conformal field theory with 
central charge c = 4/5. For N > 5, there is a critical 
"Villain" line. This line appears by arguments famil- 
iar from the Coulomb gas approach to the Kosterlitz- 
Thouless transition in the XY model. Namely, one as- 
sumes that the discrete clock variable a(x,y) can be re- 
placed by a bosonic field <p(x, y) taking continuous val- 
ues: the effect of the original discreteness is a potential 
cos(N(p) (the operator {N, 0} in the notation of section 
2). When N > 5 there is a regime in which both vortices 
(operators {0, nj for n integer) and the cos(iV</>) potential 
are irrelevant |2(J . In this regime the free boson fixed-line 
is therefore stable. In the full parameter space of the gen- 
eralized clock model with no chirality, there is a region 
which flows to this stable fixed line. However, there is 
an additional fixed point in this parameter space, known 
as the Fateev-Zamolodchikov or parafermion fixed point 
pl|p2| . This unstable fixed point is at the edge of the 
region of stability of the Villain line, and is at a self- 
dual value of the couplings. The parafermion conformal 
field theory describing this fixed point has central charge 
c = 2(N - 1)/(N + 2) while the Villain line has cen- 
tral charge c = 1. As N gets higher, the structure gets 
more elaborate, but the parafermion and Villain points 
are the only ones we need to discuss here; see |^3j for a 
thorough discussion of the phase structure of the Z5 and 
Z 6 models. 

The case of interest for the coupled Luttinger liquids 
turns out to be N = 4. The existence of the Villain 
line follows directly from the fact that the (non-chiral) 
Z4 clock model is equivalent to the Ashkin- Teller model. 
The Ashkin- Teller model is two coupled Ising models, 
with s, a = ±1 on every site related to a via 

ina/2 _ S + i<J 
l + i 

so at A = A = the action becomes 

S = — } j (Kl(<Ti<Tj + S l S J ) + K2(Ti(JjSiSj) . (25) 

<ij> 

The original chiral clock model is at K 2 = 0, where 
the two Ising models decouple. On the self-dual line 
smh(2-Ki) = exp(— 2K 2 ), the Ashkin- Teller model is, 
in turn, equivalent to a special case of the lattice 8- 
vertex model which is not only integrable but critical 
for K 2 < ln(3)/4 Q. The continuum field theory de- 
scription of this critical line is precisely a free boson (or 
to be precise, its orbifold pq|). That the Ashkin- Teller 
model is equivalent to a free boson should not be sur- 
prising. The two-dimensional Ising model is equivalent 
to a free Majorana (real) fermion, so two Ising models are 
equivalent to a Dirac fermion. In the continuum theory, 



the four-fermion coupling is the only way to add inter- 
actions without moving off the critical line. The N — 4 
parafermion critical point is a special point on this line, 
namely K 2 = (ln(\/2 + l))/4. The critical point of the 
four-state Potts model (which has a full S4 symmetry) is 
on the self-dual line at K 2 = K\ = (ln3)/4. 

Now that we have found a critical line of the lattice 
model and the corresponding field theory, we need to 
find the operators in the field theory which are relevant 
perturbations. The thermal perturbation is not chiral: it 
corresponds to moving off the self-dual line while keeping 
A = A = 0. This model is not integrable on the lat- 
tice, but the dimension of the thermal operator is known 
by combining the exact solution at T = T c 0] with 
standard scaling assumptions and a Coulomb-gas analy- 
sis |2(|. We can put this result in our notation for free 
bosons by defining the radius tat via 

The thermal operator then is {0,±2}at, with left and 
right dimensions (2r\ T , 2r 2 AT ). At the K 2 = decoupled 
point, tat = 1/2, so the thermal perturbation has di- 
mensions (1/2,1/2): it corresponds to giving a mass to 
the Ising fermions. At the parafermion fixed point, the 
thermal operator is of dimension (1/3, 1/3). 

We find the operators in the long- wavelength field the- 
ory which couple to A and A by following the argument 
in JIt]]. When A and A are small, the linear term in (|24| ) 
can be rewritten as 

K X (Asm[^f(D x a)}+A S in[^i(D y a)]j . 

As one would expect from difference operators, smD x f 
and sin D y f transform into each other under 90° lattice 
rotations. In other words, these operators are not invari- 
ant under Lorentz transformations: they transform like 
vectors. In the continuum theory, they should therefore 
have Lorentz spin ±1 + in, for some integer n. The chi- 
ral perturbation must correspond to adding to the action 
operators with Lorentz spin ±1 (operators with higher 
spin are irrelevant). We can identify these operators un- 
ambiguously by using one other fact about this pertur- 
bation: second-order contributions of the chiral terms 
renormalize the original thermal perturbation, because 
2 sin 2 x = 1 + cos 2x. Thus the operator product of the 
chiral operator (Lorentz spin I) with the antichiral oper- 
ator (Lorentz spin —1) must contain the thermal opera- 
tor. Under operator products, the electric and magnetic 
charges merely add. This means that the chiral operators 
must be {1, 1}at and { — 1,-1}^, while the antichiral 
operators are {1, — 1}at an d {lj — 1}at- 

For example, at the point K 2 = where the Ising 
models decouple, the interaction between Ising spins con- 
nected by a bond in the x direction is 

Ki cos A(erjO- l+ i + SiS i+ i) + K\ sin A(a l s l+ i - s l a l+ i). 



G 



The second term is thus proportional to aD x s — sD x a. 
The continuum field theory therefore has action 

S=Sp+A L J (a(x,t)d L s(x,t) - s(x,t)d L a(x,t)) 
+ A R J (a(x, t)d R s(x, t) - s(x, t)d R a(x, t)) 

where Sj2 is the action of two decoupled Ising models, 
and Ax,x{ = (A ± iA). We have continued to real time 
t via y = it. The operator product of two magnetization 
operators in the Ising model indeed contains the ther- 
mal operator, so the second-order contribution of this 
term indeed obeys our above condition. Since the mag- 
netization fields cr, s of the Ising models have dimensions 
(1/16, 1/16), the operator coupling to Ax, has dimensions 
(9/8, 1/8), while the operator coupling to A R has dimen- 
sions (1/8,9/8). These operators are not total deriva- 
tives, so in the bosonic language they do indeed corre- 
spond to the operators {±1,±1} at and {±l,=Fl}^j' at 
the decoupling radius tat = 1/2. 

The action of the continuum field theory describing the 
4-state chiral clock model near its critical line is therefore 
found by adding all the different perturbations to the free 
boson action Sq: 

S = S + (T-T C )J ({0,2} + {0,-2}) 

+Ax J ({1, 1} + {-1, -1}) + A R J ({1, -1} + {-1, 1}). 

All the operators are defined with radius tat ■ This is the 
action of our coupled Luttinger liquid written in terms of 
9?_ ! One needs to make the identification t^ R = Ax,x{, 
T — T c = V%kp i and tat = \/s/2. The critical line 



disordered 



and hence the description is valid for 



„4if 2 



< 3, which 



corresponds to 1/8 < r AT < 3/8. Thus the coupled 
Luttinger liquid is equivalent to the continuum limit of 
the generalized chiral clock model if 1/4 < g < 3/4. 
For the dual regime 4 > g > 4/3, one instead identifies 
tat = VV^g and T-T c = J. 



V. COUPLED LUTTINGER LIQUID PHYSICS 
FROM THE CHIRAL CLOCK MODEL 



The chiral clock models were the subject of quite a bit 
of study a few decades ago; a review can be found in [15|. 
The four-state model with K-i ^ was discussed in [27|. 
When A = 0, the four-state model is believed to have a 
phase diagram which looks like Fig. 1. 




FIG. 1. The phase diagram of the four-state chiral clock 
model at a fixed value of K^jK\ and A = 0; Ai is given by 
(pti|). The field theory description should be valid in some 
region around T — T c , A = 0. 

This picture is for some fixed value of K<xjK\ (fixed g 
in the Luttinger problem). The critical point at A = 
is the corresponding critical point on the Ashkin- Teller 
self-dual line; changing the temperature corresponds to 
changing and K\ while keeping K%jK\ constant. 
This phase diagram should be valid for 2 — \/3 < g < 3/4. 
Since everything is invariant under the transformations 
g — ► 1/g, the models are also equivalent for 4/3 < g < 
2 + V3- For g < 2 - y/3, the chiral perturbation is irrel- 
evant, so there is no incommensurate phase. The other 
bound of the map, g — 3/4, is the value Ki = — oo. 
As shown in section 3, there is definitely an incommen- 
surate phase at the free fermion point g = 1. As we 
will argue below, we expect three types of behavior. For 
g < 2 — ^/i and its dual, the model is always commen- 
surate. For 2 — \/3 < g < 3/4 and its dual, both phases 
are possible, depending on the relative strength of the 
interchain interactions and the interchain hopping. Be- 
cause as one brings g closer to 1 the hopping operator be- 
comes more relevant and the non-chiral operator less rel- 
evant, our guess is that the incommensurate phase takes 
a larger portion of the phase diagram as g — > 3/4. For 
4/3 > g > 3/4, we argue that the commensurate phase 
vanishes. These possibilities are summarized in figure 2. 
Explaining these phase diagrams is the purpose of this 
section. 
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FIG. 2. The phases of the coupled Luttinger liquids as a 
function of the intrachain interaction strength g. 

The evidence for the existence of an incommensurate 
phase in the chiral clock model derives from several argu- 
ments. In the parity-invariant model (A = 0), the main 
evidence comes from looking at the lattice model close 
to zero temperature (K\ — > oo, K2 — > oo with K2/K1 
fixed). At zero temperature, the system is completely 
ordered, except at special points at which there are two 
lowest-energy states. In the four-state model, this hap- 
pens at A = A/ < 1/2, where 



— = V2sm 



T 1 A 1 
- ( A/ - - 

2 V 2 



(26) 



At zero temperature, if A < A/, a takes the same value 
at every site; but if Aj < A < 2 — A/, D x a = 1, i.e. 
a increases by one (cyclically) every site as one moves 
to the right. These are commensurate phases. At A = 
A/, it is equally favorable energetically for D x a to be 
or 1. Thus there are a huge number of ground states, 
characteristic of an incommensurate phase. At A = A/ 
(while still A = and T = 0), the system forms stripes 
with constant a in the y-direction. |l4| . The values of a 
increases (cyclically) from stripe to stripe as x increases. 

We can extend the argument of Jl4|] to K 2 ^= to 
show that there is an incommensurate phase at small 
but nonzero values of temperature. At non-zero tem- 
perature, the stripe boundaries are no longer necessarily 
strictly parallel to the y-axis: they can develop kinks. 
The partition function can be computed at low tempera- 
tures if one neglects stripe creation (i.e. stripes which do 
not exist for all values of y) . Defining p to be the number 
of stripes per unit length, one finds 



(27) 



When the magnitude of the right-hand side is greater 
than one, the phase is ordered and commensurate. In 
the range of K\ and A where this magnitude is less than 
1, the density of stripes varies continuously with the pa- 
rameters. This is the incommensurate phase. 

The low-temperature approximation thus accounts for 
the bottom part of the phase diagram. We now consider 
the evidence supporting the upper part of the phase di- 
agram of figure 1. First of all, it is useful to consider 
what happens for clock models for N > 4. As noted 
before, there is a region of the phase diagram in which 
vortices (which correspond in the low-temperature pic- 
ture to stripe-creation operators) and the cos N<p clock 
operator are irrelevant. The chiral operators coupling to 
A still remain relevant. Thus even at A = there is 
a range of temperatures with a phase intermediate be- 
tween commensurate and disordered, in the manner of 
p8j. It is natural to assume that as one turns on A, this 
intermediate phase turns into the incommensurate phase. 
Thus the phase diagram for N > 4 looks like figure 1 , ex- 
cept that the incommensurate/disordered transition line 
hits the y-axis above the commensurate/incommensurate 
line. For N = 4, there is a particular point on the 
y-axis at which the vortex and the clock operators be- 
come (naively) marginal, while, as we have seen, the chi- 
ral operators remain relevant. This results in the pic- 
ture of figure 1, in which the disorder/incommensurate 
and commensurate/incommensurate lines hit the y-axis 
at the same point. The system can be tuned across the 
disorder/incommensurate boundary by varying the chiral 
and non-chiral perturbations, i.e. by changing the inter- 
chain tunneling and interaction strengths in the coupled 
Luttinger liquid model. 

As g increases to 3/4, the size of the incommensurate 
region increases because the exponent in ^% decreases. 
In fact, as g — + 3/4, K% — ► —00 while K2/K1 — > — 1. 
This means effectively that the critical point at T = T c 
is moving towards the T = axis. In addition, the in- 
commensurate value A/ also is moving to 0. Moreover, 
at g = 3/4, 2(i^i + K2) = In 2 on the self-dual line so 
the exponential factor in ( p7|) goes away. Thus the com- 
mensurate phase in the field theory is being pushed out. 
Therefore, we believe that in the low-temperature phase 
at g = 3/4 the model is gapless even at A = (and 
that the low-energy fixed point turns out to be the three- 
state Potts critical point). There is no longer a com- 
mensurate phase in the field theory. We believe that for 
4/3 > g > 3/4, only the incommensurate and disordered 
phases are possible, which dovetails with our analysis of 
the weak-coupling RG equations (|2^). These equations 
are of Kosterlitz-Thouless type, and there is no evidence 
for any commensurate phase. 
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Thus we see that the g « 1 weak-coupling picture 
should persist until g is lowered to 3/4. For 2 — \/3 < 
g < 3/4, the phase diagram in Fig. 1 is applicable. In the 
continuum picture, this phase diagram results from the 
interplay between the chiral and non-chiral operators. At 
T = T c and A = 0, the coefficient of the thermal (inter- 
chain interaction) operator is tuned to be zero, and we 
have a critical point. Taking A > while keeping |V2& F | 
small then results in the incommensurate phase. Thus 
we believe that the phase structure is determined by the 
relative strengths of the tunneling and the interactions. 
(We will present additional evidence for this below.) We 
expect the size of the incommensurate phase to get bigger 
as g increases to 3/4, because of the low-temperature ar- 
gument above, and because the chiral tunneling operator 
is getting more relevant while the non-chiral interactions 
are getting less relevant. 

The incommensurate phase here is believed to be gap- 
less. The main evidence for this comes from a closely- 
related model: the sine-Gordon field theory in a back- 
ground field coupling to the charge [ p9|j30| ]. This model 
clearly has a phase diagram similar to that of figure 1. 
The commensurate phase comes when the cos( Ncj>) term 
is relevant and the background field is small. For large 
enough background field, it is energetically favorable to 
have kinks in the ground state, even though they are 
gapped. This results in the incommensurate phase, which 
clearly is gapless, because increasing the background field 
continuously increases the number of kinks in the ground 
state. 

Note that the disordered phase can have minimum en- 
ergy excitations which are commensurate or incommen- 
surate. The former case occurs in the limit in which the 
non-chiral perturbation is much larger than the chiral 
perturbation. The latter was discussed in the context 
of the weak-coupling analysis at the end of section III. 
Both of these possibilities lie within the disordered phase 
of the chiral clock model and arc, hence, adiabatically 
connectable. 

Note also that our results predict the existence of a Lif- 
shitz point (i.e. a point at which the commensurate, in- 
commensurate, and disordered phases all coexist) in the 
chiral clock model for 1/4 < K 2 /K 1 < 2 - y/3 = .268 . . .. 
In this regime, the chiral operator is irrelevant near the 
A = critical point, so the disordered/commensurate 
transition line should extend out for finite A. However, 
for g > 1/4 there is still an incommensurate phase at low- 
enough temperature. The phase diagram should look al- 
most like figure 1 , except there is a direct phase transition 
between the disordered and commensurate phases for A 
small. The phase boundary between the disordered and 
commensurate phases ends at a Lifshitz point, beyond 
which these two phases are separated by the incommen- 
surate phase. 

Strong support for the above picture comes from an 
integrable model, usually called the chiral Potts model 



Jl6| . This model comprises a special subspace of the cou- 
plings in the general clock Hamiltonian (24). The name 
is somewhat misleading: the integrable model is not Sn 
symmetric but only Zjv symmetric. When there is no chi- 
ral interaction, the integrable model reduces to the Tin 
Fateev-Zamolodchikov parafermion critical point (there 
is no Potts critical point for N > 4). 

Unfortunately, except at the critical point the cou- 
plings in the integrable model do not include the value 
A = in which we are most interested. However, it does 
include a very interesting special case, corresponding to 
A = ±iA. Because of the i, the Boltzmann weights 
are not real, but the Hamiltonian in real time is Hermi- 
tian. In the language of the coupled Luttinger model, 
this corresponds to setting the tunneling at one of the 
Fermi points, say A# =t± = 0. With this condition, the 
model remains self-dual even away from the critical line 
f[l|| . The thermal operator is not self-dual, so it cannot 
be in the continuum action either |l7| . Thus the field 
theory of the integrable model on this line is equivalent 
to that of Luttinger liquids at g = 1/3 (the parafermion 
value) coupled by a single chiral tunneling operator |L7| . 
This allows us to isolate the physics of the chiral opera- 
tor (with coefficient A^ = t±) from that of the thermal 
(p2fe P ) perturbation. 

At a special value of Al known as the superintegrable 
point, the exact excitation energies can be computed 
p2[ . One finds that level crossing has occurred: there 
are states with Zm charge (not to be confused with the 
electric charge of our Luttinger liquids) which have en- 
ergy below that of the original ground state, which is 
chargeless. Thus the ground state at non-zero tempera- 
ture contains a mixture of these Zn charged states, just 
as in the sine-Gordon model at large enough background 
fields. Thus, at the superintegrable point on the self-dual 
line, the model is gapless and incommensurate. More- 
over, this phase persists even when a small non-zero An 
is allowed. 

One very interesting exact result is that the expo- 
nents are the same at the superintegrable point and the 
parafermion critical point, once the anisotropy is prop- 
erly accounted for. This is a strong indication that the 
model is gapless and incommensurate as one tunes A 
from to the superintegrable point, and in fact all along 
the self-dual line. All other exact results support this 
assertion as well [p2[ . The continuum language also 
strongly supports this picture |l7| . Because the critical 
point is perturbed only by a single chiral operator, the 
thermal operator is not induced in perturbation theory. 
The only effect of renormalization is to scale the coeffi- 
cient of the chiral operator, and not create a gap. This 
is also implied by the analysis of the sine- Gordon model 
in a background field ]30|| , which indicates that there is 
no Lifshitz point in the three-state model. 

Thus the integrable model provides strong evidence 
that adding a single chiral operator to the action results 
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in a gapless incommensurate phase. This is in perfect 
harmony with the picture that adding chiral operators 
to the action tends to cause an incommensurate phase, 
while adding non-chiral ones causes commensurability. A 
transition occurs as the relative couplings are varied. Un- 
fortunately, we do not know how to determine the exact 
location of the transition line (except at low tempera- 
tures in the lattice model). Nevertheless, we believe that 
all available evidence agrees with the phase diagram of 
figure 1. 

VI. CONCLUSIONS 



clock model). These conclusions are consistent with ear- 
lier analytical [||,0,|]|^|],[ll],|| and numerical |Q studies. 
They contradict some of the claims of although we 
do find that it is possible for c-axis transport to be sup- 
pressed in a regime in which the tunneling is still relevant. 

These results suggest the following lessons for the in- 
tcrlayer physics of the high-T c superconductors. In-plane 
non- Fermi liquid behavior alone is not likely to be enough 
to frustrate single-electron hopping between copper-oxide 
planes in high-T c superconductors. However, the forma- 
tion of a gap (possibly with interlayer correlations) could 
have this effect, thereby leading to a state which is metal- 
lic within the planes but does not show c-axis dispersion. 



Our main conclusions are 

(1) For g = 1/3, interchain tunneling alone (which can be 
disentangled from interactions by taking t\ = V2k F = 
but t± 7^ 0) leads to the formation of symmetric and 
antisymmetric bands with displaced Fermi surfaces - i.e. 
c-axis dispersion and metallic interchain transport. 

(2) Interchain interactions which couple the 2kp density 
oscillations of the two chains suppress c-axis dispersion 
and transport. As a result, for 3/4 > g > 2 — V3 (and the 
mathematically equivalent dual case 4/3 < g < 2 + ^3), 
the relative strengths of interchain tunneling and in- 
terchain interactions determine the phase boundary be- 
tween a state with c-axis dispersion and metallic c-axis 
transport and a state with a partially gapped spectrum 
(the total charge mode remains gapless) which exhibits 
insulating c-axis transport. As interchain interactions are 
increased within the gapped phase, the c-axis dispersion 
smoothly disappears; no phase boundary is crossed as 
the minimum energy excitations move to the same point 
in momentum space. 

(3) There is no fully gapless phase without c-axis disper- 
sion and metallic interchain transport, contrary to the 
claim of j| . 

These results were drawn by mapping the problem of 
coupled Luttinger liquids onto the 4-state chiral clock 
model. The phase diagram of the latter model can be ob- 
tained from an analysis of the low- and high-temperature 
expansions of this model, together with an analysis of 
the Z4 parafermion critical point of the non-chiral lim- 
iting case of the model. Result (2) follows from this 
phase diagram. Note, in particular, that this phase di- 
agram has a single gapless phase: the incommensurate 
phase. Hence, it is not possible to have a completely 
gapless phase with two identical Fermi surfaces unless 
g < 2 — v^, in which case interchain tunneling is irrel- 
evant (this case is beyond the regime of validity of the 
mapping to the 4-state chiral clock model). Result (1) 
was obtained from exact solutions of the 4-state chiral 
Potts model (a special case of the general 4-state chiral 
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